The square- lattice F model revisited: a loop-cluster 
update scaling study 

M Weigelj and W JankeJ 

fDepartment of Physics, University of Waterloo, 200 University Av W, 
Waterloo, Ontario, N2L 3G1, Canada 

JInstitut fiir Tlieoretisclie Pliysik, Universitat Leipzig, Augustusplatz 10/11, 
04109 Leipzig, Germany 

E-mail: weigeKSboromir .uwaterloo . ca, jankeaitp.uni-leipzig.de 

Abstract. The six-vertex F model on the square lattice constitutes the unique 
example of an exactly solved model exhibiting an infinite-order phase transition of 
the Kosterlitz-Thouless type. As one of the few non-trivial exactly solved models, 
it provides a welcome gauge for new numerical simulation methods and scaling 
techniques. In view of the notorious problems of clearly resolving the Kosterlitz- 
Thouless scenario in the two-dimensional XY model numerically, the F model in 
particular constitutes an instructive reference case for the simulational description 
of this type of phase transition. We present a loop-cluster update Monte Carlo 
study of the square-lattice F model, with a focus on the properties not exactly 
known such as the polarizability or the scaling dimension in the critical phase. 
For the analysis of the simulation data, finite-size scaling is explicitly derived 
from the exact solution and plausible assumptions. Guided by the available exact 
results, the careful inclusion of correction terms in the scaling formulae allows for 
a reliable determination of the asymptotic behaviour. 
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1. Introduction 

An ice-type or vertex model was first proposed by Pauling as a model for (type 
I) water ice. It was known that ice forms a hydrogen-bonded crystal, i.e., the 
oxygen atoms are located on a four-valent lattice and the bonding is mediated by 
one hydrogen atom per bond. Pauling proposed that there be some non-periodicity in 
the arrangement of the hydrogen bonds in that the hydrogen atoms could be located 
nearer to one or the other end of the bond. This positioning should satisfy the ice 
rule, stating that always two of the bonds are in the "close" position and two are 
in the "remote" position with respect to the considered oxygen atom. Denoting the 
position of the hydrogen atom by a decoration of the bond with an arrow pointing 
to the closer oxygen, this leads to the arrow configurations depicted in figure ^ when 
for simplicity placing the oxygens on a square lattice instead of the physically realized 
diamond lattice. Generalizing the resulting six-vertex model for square ice, one assigns 
energies Ci, i = 1, . . . , 6, to the vertex configurations depicted in figured resulting in 
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Figure 1. Allowed arrow configurations for tlie six- vertex model on the square 
lattice, restricted by the ice rule. 



Boltzmann factors uji = exp(— /3ei), where (3 = l/ksT is the inverse temperature 
or couphng. Assuming an overaU arrow reversal symmetry (corresponding to the 
absence of an external electric field), one abbreviates a = lui = UJ2, b = = 
and c = W5 = We- Then, the original ice model corresponds to the choice et = 0, 
i = 1, . . . , 6, whereas another especially symmetric version assumes 

Ca = Efc = 1, fc = resp. a = b = e^'^, c = 1, (1) 

which is known as the F model of anti-ferroelectrics (2], since due to the choice of 
weights the vertex configurations 5 and 6 will dominate for low temperatures, resulting 
in a ground-state of staggered, anti- ferroelectric order as depicted in figure 13 

The six-vertex model as well as the more general eight- vertex models, obtained by 
including sink and source vertices with all four arrows pointing in and out, respectively, 
have been exactly solved in zero field using transfer matrix techniques, see reference j^l- 
They exhibit rich phase diagrams featuring first-order and continuous phase transitions 
as well as multi-critical points. In particular, the six-vertex F model undergoes an 
infinite-order transition of the Berezinskii-Kosterlitz-Thouless (BKT) type to an anti- 
ferroelectrically ordered phase and the scaling behaviour of the basic thermodynamic 
quantities can be extracted from the closed-form solution. Since there is no solution 
of the model in a (staggered) field, however, information about properties related to 
the polarization is incomplete. The same is true for the correlation function, which 
can only be evaluated at the so-called free-fermion point of the model 0, 0] (the 
correlation length, however, is exactly known for all temperatures, see below). Also, 
since the solution was obtained in the thermodynamic limit, information about finite- 
size scaling (FSS) is not exact, but must be deduced from scaling arguments. Apart 
from its prominent position as a non-trivial solvable model of statistical mechanics, 
the F model has enjoyed sustained interest due to its equivalence to the BCSOS 
surface model |3, and hence several dynamical generalizations of the six- vertex 
model have been considered 0. A six- vertex model with so-called domain- wall 
boundary conditions has recently attracted considerable interest and found numerous 
applications in counting problems, the quantum inverse scattering method etc. 0. 

The Berezinskii-Kostcrlitz-Thouless scenario of an infinite-order phase 
transition induced by the unbinding of vortex pairs in the two-dimensional XY model 
has been found exceptionally hard to confirm numerically |^ 0, Ull IT^ . This is 
partially due to the nature of the infinite-order phase transition itself, which is not 
easy to distinguish from a finite-order transition numerically, and the presence of a 
critical phase, which render many of the standard FSS techniques less useful. The main 
trouble, however, is caused by the presence of logarithmic corrections, expected to be 
present on general grounds for a theory with central charge c = 1 13] and explicitly 
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Figure 2. Cutout of one of the two anti-ferroelectrically ordered ground 
states of the square-lattice F model. The state consists of vertices 5 and 6 at 
equal proportions. The dashed lines indicate one of two tilted sub-lattices of 
ferroelectrical order. 



found from the BKT theory of the model [lj| . While a numerical confirmation of the 
leading scaling behaviour of the BKT transition has been achieved in the past decade 
or so, the resolution of the logarithmic scaling corrections is still at the forefront of 
problems amenable to numerical investigation today ll^ Il7| | . 

From duality arguments and mapping to Coulomb gas systems, the F model is 
known to be asymptotically equivalent to the two-dimensional XY model at criticality. 
Thus, apart from being an interesting subject in its own right, a detailed analysis of 
the thermal and FSS properties of the six-vertex F model in the critical phase and at 
its BKT point serves as a guideline for simulations of the XY model case. Guidance 
is been given here through the fact that the exact solution of the F model yields the 
leading singularities including the correction terms explicitly, and, most notably for 
numerical purposes, the exact critical coupling of the model. Uncertainties occurring 
in analyses of the XY model such as systematic errors in the determination of the 
transition point or the effect of neglected higher-order correction terms can be studied 
rather explicitly for the F model. Finally, when it is found here that one has to 
consider large system sizes and proceed carefully when including correction terms into 
the fits, this situation should also be put into relation with the case of an F model 
placed on an annealed ensemble of random lattices considered recently [l^ . Guided 
by the present investigation, this case has to be analyzed even more carefully due to 
an additional fractality of the lattices, which reduces the effective linear extent of the 
amenable lattice sizes, thus increasing finite-size effects even further. 

The other paradigm example of an exactly solved non-trivial model of statistical 
mechanics, the two-dimensional Ising model, has served as a benchmark system and 
playground for new ideas in the theory of critical phenomena as well as for new 
algorithms in computer simulations in an overwhelming number of studies, and almost 
all of its aspects have been investigated (but not necessarily understood). In contrast, 
for the case of vertex models only rather recently efficient cluster-update Monte Carlo 
algorithms have been developed [19|, |20|, |2l| , mainly with the mapping of vertex models 
on quantum chains in mind, and some simulations of special aspects of the six- vertex 
model, such as dynamical critical exponents of the considered algorithms 20, 2^, 
properties of the equivalent surface models 23], matching of rcnormalization-group 
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flows with the XY model j2^], or the case of domain- wall boundary conditions |2l| 
have been analyzed. A systematic thermal and FSS study of the F model in the 
critical high-temperature phase, at the critical point and its low-temperature vicinity 
including the analysis of the logarithmic correction terms, however, is to our best 
knowledge lacking so far. 

The rest of the paper is organized as follows. In section 2 we outline the extent 
of exact knowledge about the phase diagram and the occurring transitions of the six- 
vertex model and the F model in particular and give an overview over scaling at a 
BKT point in general. After a short description of the simulational setup used, section 
3 contains a report of the analysis of the simulation data, comprising the FSS analyses 
of the critical-point thermodynamic properties (where the corresponding FSS relations 
are explicitly derived from the closed- form solution), an investigation of the behaviour 
in the critical phase as well as a thermal scaling analysis in the low-temperature phase 
of the model. Finally, section 4 contains our conclusions. 



2. Analytical Results 

2.1. Exact solution and phase diagram 

The square-lattice, zero-field six-vertex model has been solved exactly in the 
thermodynamic limit by means of the Bethe ansatz by Lieb j25l | and Sutherland |2^ . 
The analytic structure of the free energy is most conveniently parametrized in terms 
of the reduced coupling 

such that the free energy takes a different analytic form depending on whether A < — 1, 
— l<A<lorA>l. This leads to a phase diagram of the model consisting of four 
distinct phases as shown in figure O The phases I and II are both characterized 
by A > 1, thus corresponding to the same analytic form of the free energy; they 
represent ferroelectrically ordered phases, the ground-states being related to each 
other through a global rotation by tt/2. In these phases, the system exhibits the 
peculiarity of sticking to the respective ground-states also for non-zero temperatures. 
The intermediate case — 1 < A < 1, corresponding to phase III, includes the infinite 
temperature point a — b = c = I and thus belongs to a disordered phase, which turns 
out to be massless, i.e., it exhibits algebraic correlations throughout. This latter effect 
can be traced back to the fact that the six-vertex model corresponds to a critical 
surface in the phase diagram of the eight-vertex model. Finally, for A < —1 one 
has c > a + b, resulting in the anti-ferroelectric order of phase IV described above 
for the F model. The parameter space of the F model is restricted to the dashed 
line connecting phases III and IV indicated in figure |3| The dotted line of figure 13 
indicates the curve A = 0, where the six-vertex model is equivalent to a system of 
free fermions and an exact solution is even possible in the presence of a staggered field 
0, 0| . The nature of the transitions between the phases I-IV can be extracted from 
the exact solution 0, [i^ . Crossing the phase boundaries I — > III and II III 
one finds discontinuities corresponding to first-order transitions. The transition III — > 
IV, on the other hand, is peculiar in that all the temperature derivatives of the free 
energy exist and vanish exponentially as the transition is approached. These are the 
properties of a BKT phase transition to be detailed below in section 




Figure 3. The phase diagram of the square-lattice, zero-field six-vertex model 
in terms of the re-scaled weights a/c and fe/c. Phase boundaries are indicated by 
solid linos. The dashed line denotes the parameter range of the F model. The 
dotted line corresponds to the free-fermion line A = 0. 



While the ferroelectrically ordered phases I and II exhibit a plain polarization, 
which can be used as an order parameter for the corresponding transition, the anti- 
ferroelectric order of phase IV is accompanied by a staggered polarization with respect 
to a sub-lattice decomposition of the square lattice. This is equivalent to a mutually 
inverse plain polarization on two tilted, square sub-lattices as indicated in figure [3 
An order parameter for the corresponding transition can be defined by introducing 
overlap variables ai for each vertex i such that 

CTj = Ui * (3) 

where denotes the anti- ferroelectric ground-state configuration depicted in figure |21 
and the product denotes the overlap given by 

4 

v*v' = J2^k{v)Akiv'), (4) 
fc=i 

where k numbers the four edges around each vertex and Ak{v) should be +1 or —1 
depending on whether the corresponding arrow of v points out of the vertex or into it. 
The thus defined spontaneous staggered polarization Pq = {(Ji)/2 = (o')/2 constitutes 
an order parameter for the antiferroelectric transition III — > IV. 

As indicated above, the F model orders anti- ferroelectrically at A = —1, 
corresponding to a critical coupling = In 2. From the exact solution 0, HE l^if . 
the model's asymptotic free energy per site in the low-temperature phase /? > /3c can 
be written as 

2 TOCOsh(mA) 

m— 1 ^ ^ 

where A = acoshfi exp(2/3) — 1]. On the high-temperature side /3 < /3c it takes a 
different analytic form and has the following integral representation, 

fhigh. ^ ^g_J_ r In ( cosh^-cos2^ \ 

^ ^ 4^1 Jo cosh{TTt/2^i) y coshi-1 J' ^' 
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where fj, = arccos[i exp(2/?) — 1]. The correlation length is given by the two equivalent 
expressions 




where x = exp(— 2A) and y = exp(— 7r^/2A) are "dual", conjugate nomes of an elliptic 
function, yielding two different representations being rapidly convergent for large A 
(first form) and small A (second form), respectively. Although the general F model 
has not been solved in a staggered electric field, as a single result the spontaneous 
staggered polarization is known exactly for all inverse temperatures /3 > (3c 22}, 

Po(A)V2 ^ J|tanh(nA) = f + 2 ^(-f )".x"^ = f^j (8) 

n— 1 n—1 ^ ^ n—1 

where again the first two forms are rapidly convergent for large A, away from criticality, 
and the third form converges fast close to the critical point A = 0. As a proper order 
parameter, the spontaneous polarization Pq vanishes identically in the critical high- 
temperature phase P < Pc- 



2.2. The Berezinskii-Kosterlitz-Thouless phase transition 

As stated, the F model undergoes an finite-order phase transition of the BKT type 
at /3c = In 2. For later reference, let us shortly bringto mind the basic features of the 
BKT scenario for the two-dimensional XY model |8j, which forms the paradigmatic 
case of an infinite-order phase transition, albeit the exact solution of the F model 
was published a couple of years earlier. As a consequence of the Mermin-Wagner- 
Hohenberg theorem |23|, the two-dimensional XY model cannot develop an ordered 
phase with a non-vanishing value of a locally defined order parameter for non- 
zero temperatures|. Nevertheless, it undergoes a finite-temperature phase transition 
resulting from the unbinding of vortex pairs superimposed on an effective spin-wave 
behaviour of the low-temperature phase. Above the critical temperature, spin-spin 
correlations decay exponentially, 

G(r)~e-'-/«(^), TiT„ (9) 

while below Tc long-range correlations are encountered, 

G(r)^r-''(^), T<re, (10) 

such that the correlation length £,{T) = oo for all T < Tc and the massless low- 
temperature phase corresponds to a critical line terminating in the critical point Tc 
1^. The critical exponent rj — r]{T) varies continuously along this critical line, with 
rjc = ^{Tc) = 1/4. Approaching the critical point Tc from above, the correlation 
length diverges exponentially instead of algebraically as for a usual continuous phase 
transition, 

e(T) ~ exp(aA''), <>0, (11) 

where t — {T — T^jTc and p — 1/2. The behaviour of further observables at the 
transition point can be conveniently expressed in terms of this singularity of the 
correlation length. In particular, the magnetic susceptibility diverges as 

X(T) ~ r/'' = TiTc. (12) 

% Note, however, that on a finite lattice, the magnetization attains a non-zero value in the low- 
temperature phase, cf. reference i29il . 
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The specific heat, on the other hand, is only very weakly singular, behaving as 
(omitting a regular background contribution) 

a - r'- (13) 

Finite-size scaling analyses of the BKT transition are hampered by the occurring 
essential singularities and the presence of a critical phase. As a consequence of the 
latter, magnetic observables such as the susceptibility do not exhibit maxima in the 
vicinity of the critical point, which otherwise could be used for an estimation of the 
transition temperature from finite systems. For the same reason, also the Binder 
parameter req uires a more careful treatment than at a standard second-order phase 
transition [29l l30|. Nevertheless, the general arguments for finite-size shifting and 
rounding remain valid, such that suitably defined pseudo-critical points T*{L) for 
systems with linear extent L scale to the critical point as [sJl 

[r*(L)-Tj/r,^(lni)-i/p, (14) 

cf. equation Hll|l . since sufficiently close to the critical point the growth of 
the correlation length becomes limited by the linear extent L of the system. 
Correspondingly, ^ can be replaced by L (neglecting corrections to scaling for the 
time being) to yield the FSS law 

x{Tc,L)^L-^l'^ = L^-^% (15) 

which for 77c — 1/4 predicts a rather strong divergence. On finite lattices, the 
specific heat is found to exhibit a smooth peak, which is however considerably shifted 
away from the critical point into the high-temperature phase and does not scale as 
the lattice size is increased js^- Thus, with the main strengths of FSS being not 
exploitable for the BKT phase transition, the focus of numerical analyses of the 
XY and related models has been on thermal scaling, see, e.g., references In 
addition, renormalization group analyses predict logarithmic corrections to the leading 
scaling behaviour as expected for a theory of central charge c = 1, which have 
been found exceptionally hard to reproduce numerically due to the presence of higher 
order corrections of comparable magnitude (for the accessible lattice sizes) |llUl2lll7| . 

From the exact solution of the square-lattice F model, equations |(S1)-(|H1), one 
extracts the asymptotic behaviour in the vicinity of the critical point (3c — In 2. 
Approaching the critical point from the low-temperature side, A | 0, the singular 
part of the free energy density (O and the correlation length {T)) behave as 

/sing(A) - 4fcBTcexp(-7rVA), 
r^(A) - 4exp(-7rV2A). 

Since A goes as A (— t)^/^ for t ^ 0, this exactly corresponds to the essential 
singularity described above for the BKT transition of the two-dimensional XY model 
with p = 1/2. The specific heat has the weakly singular contribution Cy ^ £_~^ as 
expected. Concerning properties related to the order parameter, the situation for the 
F model is somewhat different from that of the XY model. The order parameter © is 
non- vanishing for finite temperatures in the ordered phase§. Thus, the corresponding 
staggered anti-ferroelectric polarizability x shows a clear peak in the vicinity of the 
critical point for finite lattices. However, in the limit L 00 the polarizability diverges 
throughout the whole critical high-temperature phase. Note that compared to the XY 

§ Note that the Mermin-Wagner-Hohenberg theorem |28j does not apply to the F model with its 
discrete symmetry. 
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model the roles of high- and low-temperature phases are exchanged in this respect, as 
expected from duality 32]. The spontaneous polarization (jHJ scales as 

Po(A) ^ A-i exp(-^V4A) ^ r'^' InC (17) 

as A I 0, implying f3/iy = 1/2. Assuming the Widom-Fisher scaling relation 
a -I- 2/3 7 = 2 to be valid||, from equations H16|l and p7ll Baxter conjectured the 
following scaling of the zero-field staggered polarizability ^7} , 

X(A) ^ exp(7rV2A) ^ a^n^', (18) 

which implies -y/i^ — 2— 77^ = 1, obviously different from the XY model result 77^ = 1/4. 
Since the whole high-temperature phase is critical, scaling of the polarizability is 
expected throughout this phase. In fact, for the free-fermion case A = or 
p = Pf = (ln2)/2, which is exactly solvable in a staggered field a logarithmic 
divergence of the polarizability is found, implying 2 — rjf = 0. More recently, the 
behaviour of rj in the critical phase of the F model has been conjectured from scattering 
methods to follow the form p^ . Is^ . Issf 

r/(A) = 7r/arccosA. (19) 

This is in agreement with the exact results for the critical F model at A = — 1 and 
the free-fermion case at A = 0. Additionally, the pure ice model at A = 1/2 is known 
to have a "dipolar" correlation function with 77 = 3 as predicted by ((T^ [s^ . Is^ . 
Note that, since the dual relation to the XY model is only valid at criticality and 
the XY model magnetization is not equivalent to the polarizability of the F model, 
this result is not simply related to the exponent 77 of the XY model in its critical 
low-temperature phase, which actually decreases as one moves into the critical phase, 
see, e.g., references [l^ls^ . 

The common occurrence of a BKT type phase transition for the XY and F 
models is no coincidence. In fact, it can be shown that the critical points of both 
models are asymptotically dual to each other [s^- This can be seen by noting that 
the Villain representation of the XY model |3^ is dually equivalent to a model of the 
solid-on-solid (SOS) type known as the discrete Gaussian model 1391 . which in turn, as 
typical for SOS models, can be mapped onto the Coulomb gas |40j. The F model, on 
the other hand, also has a height-model representation known as the BCSOS (body- 
centered SOS) model 0, which is itself asymptotically equivalent to the Coulomb gas. 
Alternatively, the stated equivalence can be seen from the loop representation of the 
0(71) vector model which for the critical 0(2) model yields a close-packed loop 
ensemble equivalent to that of the loop representation of the critical F model |42| . 
The apparent discrepancy regarding the magnetic exponents (3/v and j/i' between 
the XY and F models, on the other hand, is not an indicator of different universality 
classes of the models, but reflects the fact that the F model staggered polarizability 
is not equivalent to the magnetic susceptibility of the XY model. 

II Although the BKT transition is characterized by essential singularities and thus the conventional 
critical exponents are meaningless, one can re-define them by considering scaling as a function of the 
correlation length ^ instead of the reduced temperature t |3j . The exponents a, (3 and 7 used here 
and in the following should be understood in this sense. The exponent p, however, has its special 
meaning defined by llli . 
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3. A Loop-Cluster Update Scaling Study 

3.1. Simulation Setup 

For an analysis of the six-vertex F model via Monte Carlo simulations, a suitable 
simulation update scheme has to be devised. Since the focus here lies on the 
investigation of the vicinity of the BKT transition and the critical phase of the model, 
all local updates will suffer from the severe critical slowing down with dynamical 
critical exponent z fa 2 expected at or close to criticality. Fortunately, a fully-Hedged 
framework of cluster algorithms has been constructed for the simulation of the six- 
and eight-vertex models, mainly motivated by their equivalence with the Trotter- 
Suzuki decomposition of quantum spin chains. Here, we apply the so-called loop- 
cluster algorithm [4^ . which operates on a representation of the vertex model by 
polygons consisting of the lattice edges and induced by a stochastic breakup of the 
lattice vertices, for details see reference ji^l- For the case of the F model at criticality, 
a reduction of critical slowing down to z = 0.71(5) has been reported |43 |. 

Simulations were performed for square lattices with periodic boundary conditions, 
measurements were taken after each multi-cluster loop-update step due to the small 
autocorrelation times observed. To enable a proper FSS analysis, for the investigation 
of the BKT point two main series of simulations were performed; one around the 
peak locations of the staggered anti-ferroelectric polarizability for sizes L = 16, 24, 
32, 46, 64, 92, 128, 182, and 256 and another at the asymptotic critical coupling 
P^=ln2 = 0.6931 . . . with additional lattice sizes of L = 364, 512, 726 and 1024. To 
examine the behaviour in the critical phase, additional series of simulations have been 
performed at fixed temperatures (3 < j3c with lattice sizes identical to those at /9c ■ Per 
simulation, after equilibration a total of between 1 x 10^ and 2 x 10^ measurements 
were taken. 

3.2. Results of the Finite-Size Scaling Analysis 

3.2.1. Non-scaling of the specific heat. The specific heat is defined by 



where Vi denotes the configuration of vertex i of the lattice. It exhibits a broad peak 
shifted away from the critical point into the low-temperature phase ^^If- The essential 
singularity predicted by equation (|13() cannot in general be resolved, since it is covered 
by the presence of non-singular background terms. Thus, the non-scaling of a broad 
specific-heat peak (together with a scaling of the susceptibility or polarizability to be 
considered below) is commonly taken as a first good indicator for a phase transition 
to be of the BKT type Indeed, this is what is found from the simulation data 
as is shown in figure 0] No scaling is visible, apart from very minor deviations for 
the smallest lattice sizes and close to criticality. All data points collapse onto a single 
curve, which is identical to the exact asymptotic behaviour of extracted from the 
free energy density of equations © and © as displayed in figure ^ for comparison. 

% R ecall that the specific heat of the 2D XY model exhibits a peak in the high-temperature phase 
llj . as expected from duality. 




(20) 



(21) 
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Figure 4. Non-scaling of the specific lieat Cv of tlie square-lattice F model. 
For clarity, simulation results from only three lattice sizes are shown. The solid 
line denotes the exact asymptotic result found from the free energy density of 
equations J3 ^^'^ El- The dashed vertical line marks the infinite- volume critical 
point p^ = ln2 = 0.6931 . . .. 



In particular, at the critical point f3c = In 2 wc find for the internal energy U = (E) 
and the specific heat for the 1024^ lattice. 



Uif3c) = 0.333335(4), 
dPc) = 0.3005(15), 

in perfect agreement with the exact results U{Pc 
0.2989 0. 



(22) 



1/3 and C^iPc) = 28(ln 2)2/45 



3.2.2. The critical coupling. For an independent determination of the critical 
coupling j3c from the simulation data, we exploit the fact that the maxima of the 
staggered polarizability for finite lattices should be shifted away from the critical point 
according to the scaling relation H14|l . From finite-lattice simulations the polarization 
is determined by breaking the symmetry explicitly, i.e., if one defines S — ct^, the 
spontaneous polarization is measured as Po = (1 5^1) a-nd the polarizability is estimated 

by 

X^[{Y.^)-m?]/L\ (23) 

The peak locations of x(/3) were determined from simulations at nearby couplin gs 
by means of the reweighting technique 0| . The phenomenological theory of FSS |3l| 
implies that the polarizability x for ^ finite lattice can be expressed as 

x(/3,L) = L^/^X[L/e(/3,oo)], (24) 

where 5(/3,oo) = £,{P,L = oo) denotes the correlation length of the infinite system 
and X is an analytic scaling function (here, we omit additional irrelevant scaling fields 
representing corrections to scaling). Now, the maxima of x(/5, L) correspond to the 
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maximum of X and thus all must occur at the same value of the argument oo) 
(provided X only has one maximum), 

^ — const, (25) 



^[/3*(L),oo] 

thus defining a series of pseudo-critical temperatures (3*{L) — I3^{L). To find the 
general form of /3*{L) in the scaling region, we need to solve the expression ^ for /3. 
Inversion of the Taylor series of f of equation in powers of y = exp(— 7r^/2A) yields 

-1 



A = In 



Jr'-^r' + o(r')Jj , (26) 

and P{X) expands around the critical point A = as 

P = \n2+^X'^^X^ + 0{X^). (27) 
To leading order in both expansions one thus has via equation (|25|l 

P*{L) =f3c + ApilnACr^ ^I3c + Af^iXn^nL)"^ , (28) 

where Ap = (7rV4V2)^ Since in the FSS region ^ w L and the magnitude of the 
correction term ^"'^ in equation H26(l is relatively suppressed by a factor of 10""* 
already for the smallest lattice size L = 16 considered here, we conclude that this type 
of correction is not important at the available level of statistical accuracy. Taking 
higher-order terms of l|77|l into account, on the other hand, leads to the corrected 
scaling form 

j3*{L) = /3c + Ap{\riAnL)-'^ + Cp{\nAnL)-'^ + 0\{\n'i,nL)-^], (29) 

with C/3 = (7r^/8\/3)^. We tested fits of this expected asymptotic form to the 
simulation data for the toy model of an analytically generated series of pseudo-critical 
points P*{L) defined by equation and the exact form of the correlation length 
(0 with K = 1 and L = 16, 24, . . . , 256 as in the simulations, while taking /3c, k and 
the amplitudes A/^ and C/j as fit parameters. Already without the correction, i.e., 
enforcing C/j = 0, the critical coupling is reasonably reproduced as /3c — 0.692; the 
presence of neglected corrections shows up, however, in a fit result k 1.3. Lifting the 
constraint on the amplitude C^, one arrives at /3c — 0.6931 and k = 1.05, indicating 
perfect agreement with the input data on the level of accuracy to be expected from 
the simulations. 

Considering real simulation data, one might actually want to replace ^(/3, oo) in 
H25|l by the finite-size expression ^(/3, L) (or, equivalently, allow the constant k to 
depend on system size), introducing additional corrections not explicitly included in 
the scaling form (|24(l +. Exact expressions for the finite-size behaviour of ^ are not 
available*, however, in view of the scaling forms l|17|l and H18(l and the experience with 
various models with logarithmic scaling corrections, it seems reasonable to make the 
following general ansatz, 

e[/3*(L),L] = KL[l + A^(ln4Ki)"«] (30) 

+ Note that for the case of models with multipHcative logarithmic corrections, the replacement 
^(/3, oo) — > £,(P, L) o n the r.h.s. of 1241 has been suggested as the proper way of describing FSS 
in the first place l47l . 

* Note, however, that exact expression are available for the finite-size correlation length of the XY 
model on the strip geometry, cf. reference |48(|. 
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Table 1. Parameters of least-squares fits of the functional form 1311 to the 
simulation estimates for the peak locations of the staggered polarizability. No 
correction terms were taken into account, i.e., Bp = Cp = were held fixed 
throughout. Q denotes the quality-of-fit parameter. 



-^min 


/3c 


^0 


Q 


16 


0.73822(48) 


0.3547(62) 


0.00 


24 


0.73270(59) 


0.4533(87) 


0.00 


32 


0.73033(74) 


0.5018(126) 


0.00 


46 


0.72635(110) 


0.5912(223) 


0.46 


64 


0.72409(172) 


0.6453(385) 


0.88 


92 


0.72322(261) 


0.6668(624) 


0.78 


128 


0.72077(463) 


0.7307(1173) 


0.79 



with some a priori unknown exponent to^. Depending on the sign of lu^, this includes 
two basic inequivalent cases, namely a leading multiplicative logarithmic correction for 
> or an additive logarithmic correction for < 0. For the Ising and generalized 
(j)'^ models at their upper critical dimension, for instance, one has multi plic ative 
logarithmic corrections, corresponding to uj^ > 0, see, e.g., references [ZtL l49l Isoj . 
On the other hand, for the two-dimensional q = A Potts model ^| as well as the 
two-dimensional XY model both of which are asymptotically related to 6- vertex 
models, only additive logarithmic corrections to the finite-size correlation length occur 
at criticality. For positive w^, the replacement ^(/3,c») —^ L) in the derivation 
of the scaling of the pseudo-critical temperatures (i*{L) from and (|77|l 

produces a correction of the form ln(ln4KL)/(ln4Ki)^, whereas for integer < the 
corrections can be expanded in a power series in l/ln4Ki. To enable linear fits, we 
finally express the scaling forms in terms of L instead of 4kL, leading to the following 
scaling descriptions: 

l3*{L)=Pc + Ap{\nL)-^ + Bp{\nL)-'' + Cfj{\nL)-\ c^^ < 0, (31) 



resp. 

(3*{L)= ^3c + Ap{\nL)- 



\n\nL 
^^^^ 



> 0. (32) 



An indirect determination of the finite-size correlation length to be discussed below 
in section 2 . 31 st ronglv hints at the presence of only additive logarithmic corrections 
in H3()|l . implying uj^ < 0, but in some cases both possibilities will be considered here 
to illustrate the fact that a numerical discrimination between similar forms of the 
corrections is not at all easily possible [note that due to the extremely slow variation 
of the log-log term it might effectively be considered constant for the range of lattice 
sizes considered, which would render the form H32|l equivalent to the ansatz H31|l ]. 

The determined peak locations of the polarizability together with an example fit 
of the functional form (|31|l with omitted corrections, i.e., for Bf^ = Cjs = 0, to the 
data in the range L — Lmin = 92 up to L = 256, are shown in figure \^ The fit 
parameters of such fits, successively omitting points from the low-L side, are compiled 
in table H The strong deviations of the data from the form with Bp = C/3 = 
corresponding to a straight line in the chosen scaling of the axes are apparent from 
figure 13 Compared to the exact transition point l3c = In 2, the estimates of Pc from 
these fits are clearly too large, dropping only very slowly as points from the small-L 
side of the list are successively omitted, cf. table One might attempt to extrapolate 




Figure 5. Peak positions of the staggered anti-ferroelectric polarizability of the 
F model from simulations as a function of lattice size. The lines show fits of the 
form 13 li to the data. The dashed line corresponds to an uncorrected fit with 
Bj3 = Cjs = 0, starting from Lmin = 92. The solid curve corresponds to a fit with 
corrections and Lmin = 32. 



these results towards imin — > oo using the scahng form (|31|l : since the individual 
data are highly correlated, however, this would introduce a strong bias. Instead, we 
directly use the higher-order logarithmic corrections in the fitting procedure. Note 
that this effect of strong scaling corrections here occurs for rather large lattices, where 
for a usual continuous phase transition without logarithmic corrections the presence of 
scaling corrections usually would not be much of an issue for the determination of the 
leading scaling behaviour. Relaxing the constraints on Bp only or on both parameters. 
Bp and C/j, we arrive at the fit results compiled in table 13 It is apparent that the 
complexity of the completely unconstrained fit type is at the verge of exceeding the 
available statistical accuracy of the data, such that competing local minima of the 
distribution exist, which result in a rather discontinuous evolution of the amplitudes 
AjS, Bfj and Cp as the lower-end cut-off Lmin is increased. This functional form fits 
the data very well, however, and the estimates for j3c are all in agreement with the 
asymptotic value /3c = 0.6931 ... in terms of the statistical errors. It seems clear 
that the remaining vague tendency of the fits to yield /3c slightly above its asymptotic 
value could in principle be removed by including further correction terms for the 
case of extremely accurate data. Adding a (lnL)~^ term in (|3H) . for instance, yields 
j3c = 0.674(45), Q = 0.72 when including all data points. 

Using the alternative fit form (|32|l valid for > 0, on the other hand, results in 
fits quite similar to those obtained from the ansatz (|31() with both correction terms 
present. This relates back to the remark concerning the extremely slow variation of 
the log-log term, which makes it plausible that the considered correction terms can 
be effectively interchanged. No specific drift is observed on increasing the cutoff imin- 
For Lmin = 64 we arrive at (3c = 0.695(40), Q = 0.92. Hence, from the scaling of the 
pseudo-critical (inverse) temperatures (3*{L), a conclusion about the sign of can 
hardly be drawn. 
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Table 2. Parameters of fits of the form 1311 to the peak locations of the 
polarizabihty. As indicated, Cp = was hold fixed for the fits shown in the 
upper part of the table, while both parameters, and Cp, were allowed to vary 
in the fits presented in the lower part. 



-^min 


/3c 


^3 






Q 


16 


0.7103(17) 


1.58(7) 


-2.91(17) 


[0.0] 


0.83 


24 


0.7118(29) 


1.50(14) 


-2.79(37) 


[0.0] 


0.78 


32 


0.7069(46) 


1.78(25) 


-3.47(67) 


[0.0] 


0.97 


46 


0.7108(85) 


1.53(51) 


-2.74(149) 


[0.0] 


0.97 


64 


0.714(16) 


1.34(103) 


-2.16(320) 


[0.0] 


0.89 


16 


0.7119(84) 


1.44(71) 


-2.2(34) 


-0.9(46) 


0.73 


24 


0.695(16) 


3.2(16) 


-11.4(848) 


12.(124) 


0.84 


32 


0.723(32) 


-0.04(343) 


6.7(191) 


-15.8(297) 


0.96 


46 


0.713(76) 


1.3(93) 


-1.1(553) 


-2.7(920) 


0.88 



3.2.3. The correlation length. Owing to tlie original relation of the present work to 
an investigation of the F model on dynamical random graphs where the definition of 
connected correlation functions is found to be highly non-trivial [l3 |. we have not 
measured the correlation length directly. It turns out, however, that an indirect 
determination is possible. Due to the scaling forms IjlTI) and H18|) of the polarization 
Pq and the polarizabihty x, for the combination x/Pq the multiplicative logarithmic 
corrections cancel such that to leading order 

xiP, L)/P^{P, L) = A^/p.aP, Lf (33) 

in the scaling region, i.e., for (3 — P*{L). This relation allows for an indirect 
determination of the exponent uj^ of the logarithmic correction of the scaling form 

X{P*,L)/P^{(3*,L) = A^/pgL'[l + M\nLr^]\ (34) 

where, again, the dependence on k has been dropped since its inclusion leads to 
very badly converging, unstable fits. This corresponds to the omission (for the time 
being) of higher-order corrections to scaling. From the simulation data, we find the 
combination (|34|l at the critical point (3c to be very well described by a quadratic 
behaviour in L, the correction in square brackets being quite small in absolute terms, 
cf. figure Fitting the form (|34|l to the data, we find clearly negative correction 
exponents which, however, strongly depend on the cut-off Lmin, systematically 
dropping from uj^ = —0.77(96) for Lmin = 16 to, e.g., uj^ = —4.1(74) for Lmin = 64 
with qualities Q > 0.8. The large statistical errors on the estimate support the 
visual impression from figure El that the correction is actually too small to be reliably 
resolved at the present level of accuracy, whereas the systematic drift results from the 
omission of higher-order corrections. In fact, if we assume the familiar power-law form 
with negative exponents only, 

x{f3*,L)/PS{P*,L)=A^/PgL'[l + A^{lnL)-'+B^{lnLr']\ (35) 

we find stable and high-quality (linear) fits for the amplitudes as L„iin is varied, cf. 
the fit with L^i^ — 16 and Q = 0.90 displayed in figure E| Thus, in analogy to 
the cases of the two-dimensional XY [l^ and q = A Potts 0| models, the F model 
finite-size correlation length appears to exhibit only additive logarithmic corrections 
corresponding to < in H30|l . 
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Figure 6. FSS of the combination x(/3i ^)/^o at tlie asymptotic critical 

point /3c = 0.6931 . . . scaled by the leading behaviour to exhibit the corrections 
to scaling. The solid line shows a fit of the form 1351 to the data. 



3.2.4- P'SS of the spontaneous polarization. To derive FSS for the spontaneous 
polarization, consider the second form given in equation ||SJ|, P^^^ — (2tt/ X)^/'^[y^/'^ + 
y^l^ + •■•], which is rapidly convergent in the scaling window. The sub-leading terms 
in y are strongly suppressed in the scaling regime and can be neglected compared to 
correction terms to follow. Using again the expansion of y in terms of ^ = ^(/3,(X)) 
and equation (|26|l for A(^), one arrives at 



Po(/3,(X3)=Ap„(4e)-^/''(ln4er-o 



-2 



4 (40-2 
3 



(36) 



3''^' '3 In 4^ 

with Ap^ — 4/7r, ^jv =1/2 and top^ — 1. Note that this (exact) form does not contain 
any corrections of the log-log type present in the XY model correlation function [l3 |. 
To test the sufficiency of this approximation, we again use an analytically generated, 
"artificial" series of scaling data, evaluating Po(/3, oo) exactly from (jS)) for the series 
of pseudo-critical temperatures defined by (|25|l for k = 1 and the exact expression (jT)) 
for ^ = ^(/3, oo). Fitting the form H36|l without the corrections in square brackets to 
these data, taking Ap^ and (iju as fit parameters (holding wp^ — 1 fixed), we arrive 
at 13 /v = 0.5001, which is clearly sufficiently close to the exact result in terms of the 
statistical accuracy to be expected from the simulation data. We thus conclude that 
the scaling corrections in square brackets of (|35|l can be neglected for our purposes. 

Further FSS corrections arise from the behaviour H30|l of the finite-size correlation 
length. For integer < as indicated by the investigation of x/Pg above, these 
corrections can be expanded in a power series. 



-l3/u 



Pa 



1 



InL 



Cpo 



we < 0, 



(37) 



(lnL)2 (lnL)3_ 

where, again, the effect of the multiplier k is being incorporated in the correction 
amplitudes. For the case of a positive wj, on the other hand, one would find a log- log 
correction to occur, i.e., one arrives at 

InlnL" 



Po(/3*,L) = Ap„L-'^/''(lnL)"^n-'^«/2 



1 + PPo 



InL 



UJf > 



(38) 
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where additionally, the exponent of the multiplicative logarithmic correction is 
"dressed" a.s ujpg — 

As for the route in the {f3*{L),L) plane taken towards criticality, the two 
principal choices are given by the determination of Po{P* = f3c,L) at the fixed 
asymptotic critical coupling /3c = In 2 or by considering Pq [/?* = (3^ (L) , L] at the 
polarizability peak locations P^{L)jjt. Asymptotically, both approaches should give 
compatible results; the strength and composition of scaling corrections, however, might 
be noticeably different. Following the first approach, we analyze the data from lattices 
of sizes L = 16, ... , 1024. Uncorrected fits of (|S7|l with Bp^ = Cp^ = Dp^ = and 
omitted multiplicative correction term, i.e., wp^ = 0, yield exponents approaching 
the expected value logarithmically slow on successively omitting data points from the 
small-L side of the list. For L = 92, . . . , 1024, for instance, we find 13 /v = 0.4658(20), 
statistically incompatible with (3/v = 1/2. With variable wp^, on the other hand, the 
leading scaling exponent can be reasonably reproduced (Lmin = 24), 

Ap^ = 2.159(35), 
l3/v = 0.4872(76), 

ojp, = 0.109(33), ^ > 

Q = 0.14, 

but the resulting exponent of the logarithmic correction is estimated in strong 
deviation from Lupg — 1. This shortcoming can only be remedied by including the 
power-series type corrections in (I37f) . Letting only Bp^ vary, we arrive at estimates 
P/v = 0.50(60), LUpg = 1.12(14), Q = 0.13, which fit the expectations very well. In 
view of the large error estimates, however, one should not be deceived by the very 
small deviation from the exact result. In fact, even with Cpg — Dp^ ~ still fixed, 
the distribution exhibits multiple local minima and the fit results heavily depend 
on the initial parameter values. Additionally letting Cpg and/or Dp^^ vary, the fits 
get very unstable and meaningful results can no longer be found. Using constraint 
fits, however, it can be clearly seen that the small value wp,, = 0.109(33) above is 
indeed an effect of neglected higher-order scaling corrections: fixing f3/i' = 1/2 as 
well as Lupg = 1 we determine the amplitudes Apg, Bpg, Cpg and Dp,, with Q = 0.26. 
Now, fitting the form ElhiL)'" to the values of the thus determined polynomial 
l + Bp^{lTiL)~^+Cp„{lT^L)-^+Dp^{lTiL)-^, we find a « 0.83. Thus, neglecting the 
scaling corrections in square brackets of (|37() clearly leads to an effective reduction of 
the exponent estimate ujp„ from its asymptotic value ujpg = 1 by about a « 0.8 — 0.9. 

Considering the spontaneous polarization at the peak positions of the 
polarizability for lattice sizes up to L = 256, the uncorrected form with cup^ = 
Bpg = Cpg — Dpg = ylcMs very small estimates for P/i/ around (3/v m 0.25 slowly 
increasing with the cutoff imin- Including the multiplicative logarithmic correction 
of equation H37|) . i.e., relaxing the constraint up^ = 0, these results can be improved, 
and, e.g., for Lmin — 92, we find the following fit parameters, 

Ap^ = 1.49(50), 
= 0.44(11), 
u;p„ = 0.71(55), 
Q = 0.57, 

with an exponent estimate for /3/i/ well compatible with the exact result f]/^ — 1/2, 
although endowed with an unpleasantly large statistical error. The inclusion of the 

fj Obviously, the first approach is only amenable in cases where /B^ is known a priori. 
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power-series type scaling corrections of (|37|l necessary for a full resolution of the 
corrections is not possible with the present data. Fixing (3/v — 1/2, = 1, however, 
these terms provide an excellent description for the present scaling corrections and a 
quality-of-fit Q = 0.72 is already attained for Lmin = 16 (including all three terms 
■Bpo; Cpo s-nd Dpg). In general, fits to the data at the maxima of the polarizability 
are found to be somewhat less stable and precise than those to the data at f3c, which 
we attribute to the smaller available system sizes here as well as additional scatter of 
the data due to the necessary reweighting. 

In addition to the self-contained scaling routes in the {(3*{L),L) plane described, 
for the exactly solved case considered here it is also possible to perform simulations 
for the analytical series (3*{L) of inverse temperatures defined by the relation (^5)1 
with the exact expression (UJ for ^, which yields a scenario somewhat in between the 
(3*{L) — Pc and P*{L) — I3^{L) cases. This artificial series of simulation data is indeed 
found to result in quite stable fits, such that at least the amplitude Bp^ of ||1T7|) can 
be left variable to yield l3/v = 0.51(47), ujp„ = 1.1(69) with Q = 0.66 and L,„i„ = 16. 
It should be noted that also fits of the form (|38|) with a log-log correction are possible 
with good quality for the scaling of the polarization, which only very slightly change 
the estimates for f3/v, but do not lift the estimate for ujp^ to the expected value 

~ 1. Thus, it would be hard to distinguish the forms l|S7|l and solely on the 
basis of the numerical polarization data. 



3.2.5. FSS of the polarizability. Since the staggered polarizability x is not known 
exactly, a systematic discussion of x as a function of the asymptotic correlation length 
^ is not possible. However, from Baxter's conjecture H18|l the leading behaviour is 
expected to be 

x(/3,oo)=A^^/''(lnO'^^ (41) 

with ^/v = 1 and lj^ = 2. Repeating the arguments presented above for the 
polarization, again assuming an integer < in equation (|30|) . one deduces the 
following FSS ansatz, 

-Bv Cv 



< 0, (42) 



InL (lnL)2 (lnL)3_ 

whereas lu^ > would result in a form including a log-log correction as in H38|l . We first 
investigate the simulation results at criticality, using data from lattices of sizes L = 
16, . . . , 1024. From fits of the leading scaling behaviour dropping the multiplicative 
and additive logarithmic correction terms, Wj^- = and ~ C^^ = = 0, to these 
data, we find reasonable fit qualities only when dropping many points from the small- 
L side of the size range. Successively increasing the cutoff imin a very slow downward 
drift of the estimates for j/i/ is observed. For Lj^in = 92, we arrive at an estimate 
^jv = 1.0754(22) with Q = 0.96, which is clearly incompatible with the exact result in 
terms of the statistical error. Letting lo^ vary while still keeping = = ~ 
fixed, stable and good-quality fits can be attained. For the range L — 24, . . . , 1024, 
we have 

= 1.581(31), 
j/iy = 1.0166(90), 

= 0.320(40), ^ ' 

Q = 0.78, 
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Figure 7. Finite-size scaling plot of the critical staggered polarizability x(/^c, L) 
for lattice sizes from L = 16 up to L = 1024. The solid line shows a fit of the 
functional form 1421 with B-^ = = D-^ = to the data. The abscissa has 
been re-scaled such as to factor out the leading scaling behaviour oc L"'^'^ with 
-y/u = 1.0166 from the fit E^ . 

in good agreement with the exact result "f/v = 1, however, again clearly missing the 
expected asymptotic value of the exponent of the multiplicative logarithmic correction, 
u-^ = 2. The scaling plot presented in figure |7| shows this last fit together with 
the simulation data, scaled such as to expose the magnitude of scaling corrections 
present. The asymptotic value = 2 could be recovered by including the correction 
amplitudes B-,^ , and D-^ . Letting only i?^ additionally vary, coy- is already increased 
to = 1-32(17) with j/iy = 1.0(48), Q = 0.78 and L^in = 16. The obviously 
necessary higher-order terms and unfortunately cannot be fitted any more 
with the present data, however. On the other hand, fixing again j/iy = 1 and u-^ = 2, 
the present corrections can be well described by the amphtudes B^, C^, D^, resulting 
in a quality-of-fit of Q = 0.84 already for Lmin — 16. We note that here even the 
inclusion of the (InL)^'^ term is probably crucial since the leading multiplicative 
logarithmic correction is already quadratic. 

Estimates of the maxima xiPxi-^)' ^] ^^'^ available for lattice sizes L = 16, . . . , 256. 
A reasonable quality fit of the uncorrected form H42(l with uj^ = and B^ — — 

= to these data can be produced starting from Lmin = 64, which yields an 
estimate j/i^ = 1.2788(58), Q = 0.23, lying even much further off the asymptotic 
result than in the case of the critical polarizability. Letting lo-^ vary while keeping 
B^ — — Dy^ — fixed, the estimate for -f/i^ is noticeably reduced to 7/;^ = 1.13(08) 
with Lu^ = 0.71(36), Q = 0.15, for Lmin = 46, and a further tendency to decrease on an 
increase of Lmin remains. Although here again, the need for higher-order correction 
terms is apparent, we find the data not precise enough for their inclusion. Thus, 
although both methods, consideration of the critical polarizability as well as scaling 
of the peak heights of x(-^)i yield equivalent results, we find corrections to scaling 
slightly more pronounced in the latter approach. This is partly explained by the fact 
that for xiPc L) larger lattice sizes could be considered. However, even restricting a 
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fit with ^ = = for xiP^L) to L < 256, we find witli j/i^ = 1.006(11) 
for iniin = 16 a considerably more precise result closer to the asymptotic value; 
additionally, as mentioned above, no further drift of 7/^ is noticeable there as Lmin 
is increased. For the extra simulation series at inverse temperatures (i*{L) resulting 
from equation (|25|) with k — 1, a fit with = = = and Lmin = 32 leads 
to -i/v = 1.075(66), uj^ = 1.21(29) with Q = 0.48. Inclusion of the B^, and 
terms destabilizes the fits too far, although consistency with ^ / v — 1, uj^ = 2 \s again 
found on fitting the amplitudes only. 



3. 2. 6. The scaling dimension in the critical phase. Due to the criticality of the high- 
temperature phase one expects scaling and, accordingly, FSS in the whole region 
/3 < /3c = ln2. The closed-form conjecture H19|) for the exponent 77 entails predictions 
for the FSS of Po{/3,L) and xiPiL) for (3 < /3c- In terms of the scaling dimension 
Xp — l3/v = 1 — 7/21/ 13] and the inverse temperature /3, equation H19|) reads 

xp{P < /3c) - I |arccos[l - i exp(2/5)]| , (44) 
which behaves close to the critical point /3c = In 2 as 

^p{P) ^\ + — (ln2 - /3)i/2 + o(in2 - /3), (45) 

such that Xp has a vertical tangent at /3c, implying an especially sensitive dependence 
of on scaling corrections there. It is worthwhile to notice that, although the 
correspondence between the XY and F models only applies to their critical points, 
an analogous square-root singularity of the exponent 77 of the XY model is found on 
entering the critical low-temperature phase there, see, e.g., reference ^J. The leading 
scaling behaviour of Po(/3, L) and x(/3, L) for /3 < /?c is hence expected to be 
Po(/3,L) = Ap„i---(/3), 

X(/3,i) = A;,L2-2x^(/3)^ 

The solid lines of figure |S| illustrate the predicted behaviour of these exponents in 
the high-temperature phase. As can be seen, the polarizability exponent 2 — 2xp{[3) 
crosses zero at the free-fermion coupling (if — ^ln2 and, consequently, x should 
be non-divergent below. As a result, the predicted singularity would be covered by 
non-singular background terms there, such that we restrict ourselves to the range 
Pf <P<Pc here. 

To test the form (|44|l we performed seven series of simulations at inverse 
temperatures (3 — 0.35, 0.40, . . . , 0.65 with the same series of system sizes L = 
16, ... , 1024 used at /3 = In 2. Fitting the expected leading scaling behaviour H46|) to 
the simulation data, many system sizes from the small-L side have to be dropped to 
reach satisfactory fit qualities and to account for the observed slow drift of the resulting 
scaling exponents on increasing imin, which was finally chosen to be Lmin = 182 in 
most cases, cf. the data in column (a) of table O As can be seen from the fit data 
presented in figure |S1 even with this precaution highly significant deviations of the fit 
results from equation (|44|l are observed, especially close to the free-fermion coupling 
Pf. Scaling corrections are assumed here to take the form found at criticality, i.e., 

Bp„ Cpo Dp„ 



PoiP,L) = Ap„L-"^(^)(lnL)'"^"(''' 
X(/3,L) = A^L^-^''^^^\lnL)'^^^^^ 



InL (lni)2 (InL) 
, , ^x 



InL (lnL)2 (InL)^ 



(47) 
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Figure 8. Finite-size scaling dimensions xp(P) of the spontaneous staggered 
polarization and 2 — 2xp{f3) of the staggered polarizability, respectively, as a 
function of inverse temperature /3 in the critical phase /3 < /3c. The symbols denote 
results from FSS fits of the functional forms I46i resp. 1471 writh B = C = D = 
to the simulation data. The solid lines correspond to the conjecture 1441 for the 
analytic form. The vertical dashed lines indicate the locations of the the free- 
fermion point 13 f = ^ ln2 and the critical point /3c = In 2, respectively. 



As for the critical polarization and polarizability, fits including all of the correction 
terms (amounting to six variable parameters) are not possible with the available data. 
Including only the multiplicative logarithmic correction with variable exponent up^ 
resp. w^, we arrive at largely improved estimates for the scaling dimension xp in 
agreement with the prediction (|44|l . see the column (b) of tableland the data in 
figure |S1 The values of the correction exponents ujpg resp. lu-^, however, again have to 
be considered as effective exponents owed to the omission of the additive corrections 
B, C, D. Note that in principle also the values of wp^ resp. uj^ could depend on the 
value of the coupling /3 as indicated in equation H47I) . To investigate this possibility, we 
performed fits with the leading scaling exponents fixed to the presumably exact values 
of xp from H44I) . letting ujp„ resp. uj^ vary and including two orders of additive scaling 
corrections, i.e., enforcing D — only. The results of these fits are collected in column 
(c) of table|21 In all cases with the exception of /3 = 0.45, which seems to be an outlier, 
we find very good fit qualities, again indicating consistency with the conjecture iBH) . 
The estimates for topg are all consistent with a constant value of wpg = 1, independent 
of the coupling /3. The estimates for uj-,^, on the other hand, are clearly larger than 
= 2, but no general trend on varying /3 is observed. This deviation of is found 
to disappear upon inclusion of the next-order correction amplitude for which, 
however, both exponents, xp and uj^, have to be kept fixed. This term has to be 
included here but not for the polarization since for x already the leading multiplicative 
logarithmic correction is quadratic. When fixing both exponents, xp and u, fits of 
good quality are attained for both observables and all couplings /3 on including all 
three additive correction terms of H47|l . In passing we note that an analysis of x/Pq 
in the critical phase yields negative values of everywhere and fits of the power-law 
form (p!5)l describe the corrections extremely well. 
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Table 3. Fit parameters for Pq and x in the critical phase Pf < P < Pc- (a) Fits 
of the form (b) Fits of the form g3 with B = C = D = 0. (c) Fits of the 

form B71 with D = and xp fixed at the values BH . (d) Fits of the form E71 
with D = and a;pp = 1 resp. u!^ = 2 fixed. 
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3.3. Results of the Thermal Scaling Analysis 

The discussed FSS of the critical polarization and polarizability is independent of the 
value of the critical exponent p. For the scaling of the polarizability peak positions 
in section 13. 2. 21 on the other hand, the need to resolve the present strong logarithmic 
scaling corrections did not allow for an additional independent determination of p. To 
directly verify the exponential type of the observed divergences and to estimate the 
parameter p, one should hence consider thermal instead of finite-size scaling. Figure 
shows an overview of the temperature dependence of the staggered polarizability 
for different lattice sizes. The clear scaling of x for the high-temperature region 
/3 < /3c = ln2 illustrates again the presence of a critical phase. In contrast, for the 
low-temperature phase to the right of the peaks, the polarizability curves essentially 
collapse and only start to disagree as the correlation length reaches the linear extent of 
the considered lattice. Therefore, a thermal scaling analysis must be performed in the 
low-temperature vicinity of the critical point, the behaviour in the high-temperature 
phase being completely governed by finite-size effects. Here, we do not consider the 
scaling of the correlation length itself, but instead analyze the thermal scaling of the 
spontaneous polarization and the polarizability for a single lattice of size L = 256. 
Simulations were performed for a closely spaced series of temperatures in the low- 
temperature vicinity of the critical point. 
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Figure 9. Scaling of the polarizability peaks from simulation data. The 
lines simply connect the data points and are drawn to guide the eye. The 
dashed vertical line indicates the location of the asymptotic critical coupling 
Pc = 0.6931 . . . 



3.3.1. Scaling of the spontaneous polarization. From the leading term of (jSl in j/ and 
the dependence of A on (3, the spontaneous polarization behaves as 



(/3 - /3c) 



-1/2 



1/2 



X exp 



8V2 



6 



Pc) 



1/2 



(48) 



as (3c is approached from above. Taking only the leading-order terms into account, we 
consider the following scaling form, 



lnPo(/3) = Ap, + Bp^ip - Pc)-P + Cp, ln(/3 - /?,), (49) 

with Cp,, = —1/2 and p = 1/2. The window of validity of (|49|) for the thermal scaling 
of Pq for a finite lattice is limited for small deviations {3 — (3c by finite-size effects 
and for large deviations /3 — /3c by the higher-order corrections to scaling indicated in 
H48|l . If correlation lengths are measured, one might monitor the effect of the finite 
lattice size by comparing the value of the correlation length ^(/3, L) at a given (3 > Pc 
with the linear extent L of the lattice E3|tt- Here, the onset of finite-size effects is 
estimated by the beginning of the rounding of the exponential decline of Pq as /3c is 
approached. From monitoring the quality-of-fit parameter and estimation of the onset 
of the finite-size rounding, we determine a fit range of /3min = 0.77 < P < 0.85 = /3max. 
We find fits of the full five-parameter family H49|l of functions to the data possible, 
but the resulting fit parameters are endowed with astronomic error estimates and the 
corresponding distribution has multiple minima such that different "solutions" can 
be found. We thus fix one or two of the parameters at their expected asymptotic 

ft Although the behaviour of the finite-size correlation length has been indirectly analyzed above in 
section 13.2.31 unfortunately we do not have access to the amplitude to find the absolute values of 
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Table 4. Parameters of fits of the form I49i for Pq (upper part) resp. the form 
l5Ui for X (lower part) to the simulation data. Values in square brackets indicate 
that the corresponding parameter was held fixed in the fit procedure. 
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values to reach more stable fits, cf. the fit data collected in table ^ Note that the 
parameters of the fully unrestricted fit where found starting from the parameters of 
one of the restricted fits, thus explicitly selecting one of the minima. Figure [TUI 
shows the simulation data together with this unrestricted fit and the exact asymptotic 
polarization of ijHl). The vertical line denotes the inverse temperature 13* where the 
asymptotic correlation length ^(/3*, oo) of equation Q reaches the linear size L — 256 
of the system. As expected, this point approximately coincides with the inverse 
temperature where the simulation data deviate from the asymptotic result due to 
finite-size effects, thus justifying the method of determining /3min- The inset of figure 
^] shows the approximation of H48|l with only the first-order terms of both expansions 
being kept in comparison to the full asymptotic result ((HJ and the simulation data. 
As can be seen, even in the scaling range considered here, the deviation is much larger 
than the statistical errors of the data. The observed shift, however, can be mostly 
reproduced by slight changes of the amplitudes Ap^ and -Bp^, such that (|49|l still fits 
the data well. The fitted amplitudes Ap^, Bp^ and Cp^ must be considered effective, 
however, and deviations of the fitted parameters from the exact asymptotic values are 
due to the effective inclusion of neglected higher-order correction terms. 

3.3.2. Scaling of the polarizability From the conjecture (|18(l for the near-critical 
polarizability, we expect x(/3) to scale analogous to the polarization, 

\nxif3)^A^ + B^i/3-P,)-P + C^\n{P-/3,), (50) 

where the differences to the scaling of Pq only show up in the amplitudes A^, 
and = — 1. From the flattening out of the exponential divergence near Pc and by 
monitoring the quality of fit, we estimate the same scaling window /3min — 0.77 < 
(3 < 0.85 — /3max for H50|l we encountered for the polarization. We find fits for the 
polarizability to be considerably less stable than those for the polarization, and we 
did not succeed to fit all five parameters independently. Fixing = — 1, a reasonable 
result for p cannot be found, even when additionally fixing /3c = In 2, cf. the data 
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Figure 10. Thermal behaviour of the spontaneous staggered polarization Pq 
close to I3c from simulations of a 256^ system. The solid line denotes the exact 
asymptotic result JSJ, the dashed line is a fit of the form 1491 to the data. The 
vertical line denotes the point where 5(/3*,oo) = L = 256. The inset shows the 
exact solution JSJ compared to the first-order approximation of I48i in the inverse 
temperature regime used for the fit. 

compiled in tabled] Since a fit with only (3c fixed yields « 0, corresponding to 
an omission of this correction term, we also tried fits with C^^ = fixed, which work 
considerably better than fits with = —1. However, still meaningful results for p and 
can only be found when fixing one of the two parameters, which then yields good 
agreement with the asymptotic result. Figure ITTI shows the simulation data together 
with a fit with = and /3c = In 2 fixed. Comparison of the asymptotic correlation 
length (O with the system size L = 256 indicates the approximate onset of finite-size 
effects as the critical point is approached. 

To see in how far it is possible to distinguish the occurring essential singularity 
from a conventional power-law behaviour, we also performed fits to the form (|50|l with 
the left side replaced by x(/3) instead of lnx(/3) and ^ held fixed. With this 
power-law form and the same range of inverse temperatures used for the exponential 
fits, we arrive at the following parameters, 

A'^ = 13.5(49), 
B'^ = 0.040(35), 

/3c = 0.7112(96), (51) 
p' = 3.54(51), 
Q 0.13, 

where p' now would correspond to the conventional critical exponent 7 for the case 
of a finite-order phase transition. Thus, in agreement with the experience from the 
two-dimensional XY model, power-law fits can be performed with satisfactory quality 
if one accepts "unnaturally" large exponents such as p' = 3.5 here. 
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Figure 11. Thermal scaling of the polarizability on a L = 256 lattice. The solid 
curve shows a fit of the functional form l5Ui to the data, where the parameters 
Cx = and /3c = In 2 were kept fixed. The vertical dashed lines indicate the 
window of data points included in the fit. To judge the onset of finite-size 
effects, the dashed-dotted curve shows the logarithm of the ratio ^(/9, oo)/L from 
equation J7J, such that strong size effects are expected to appear as soon as 
ln[5(/3, oo)/Ll > (right scale). 



4. Conclusions 

We have considered the behaviour of the six- vertex F model on the square-lattice at its 
Berezinskii-Kosterlitz-Thouless (BKT) point and within the critical high-temperature 
phase with a series of cluster-update Monte Carlo simulations and subsequent finite- 
size and thermal scaling analyses. Due to the presence of strong logarithmic corrections 
indicated by the exact solution and expected for a theory with central charge c — 1, 
the scaling analysis has to carefully take correction terms into account and/or treat 
the presence of (even higher-order) corrections by omission of simulation points close 
to the border of the scaling region. Although the usefulness of the finite-size scaling 
(FSS) technique has been called into question at a BKT point due to the occurrence of 
essential singularities and most studies of the XY model case solely consider thermal 
scaling instead we find a FSS analysis for the F model well possible and useful, as 
long as corrections to scaling are thoroughly included. The full FSS forms including 
the correction terms are explicitly derived from the exact results augmented by the 
plausible assumption (|3U|I about the scaling of the finite-size correlation length. The 
latter is being confirmed by the analysis of a combination of observables proportional 
to a power of the correlation length without multiplicative logarithmic corrections, 
providing evidence that the finite-size correlation length exhibits additive logarithmic 
corrections in the present case (as opposed to multiplicative logarithmic corrections 
such as, e.g., at the upper critical dimension 01 )■ -Due to the ambitious nature of 
many of the fits involved, however, one has to cope with the occurrence of competing 
local minima of the distribution and a distinctive flatness of these minima in some 
parameter directions entailed by the slow variation of the logarithmic terms. We would 
like to stress that the quality-of-fit parameter Q is found to be not always sufficient for 
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the detection of neglected higher-order corrections. Omitting the discussed correction 
terms, however, the resulting estimates do not even satisfy moderate expectations 
of accuracy and are strongly biased. For the FSS analysis the knowledge of the 
exact asymptotic critical coupling /3c turns out to be highly beneficial and the results 
found from the scaling at effective pseudo-critical points are much less accurate. This 
might be taken as a caveat for simulations of the XY model, where /3c is not exactly 
known. The correction exponents upg for the polarization and uj^ for the polarizability 
could not be consistently and accurately determined in fully unrestricted fits, although 
constrained fits including further correction terms allow to establish consistency with 
the analytical solution. This experience is shared with simulational studies of the XY 
model 0, 01 . A thermal scaling analysis of the low- temperature approach towards 
criticality does only lead to reasonably precise results for the present data if at least 
one of the fit parameters is fixed to its exact value. A conventional algebraic singularity 
also fits to the data, but only when unusually large exponents are accepted. 

In addition to the analysis at criticality, we consider the scaling of the polarization 
and the polarizability within the critical high- tem pera ture phase. We find overall good 
agreement of the outcome with a conjecture [s^ l35l| for the behaviour of the scaling 
dimension xp{(5) of the polarization in the critical phase, although the resolution 
of scaling corrections appears to be even more involved here than at criticality. 
Close to the critical point scaling corrections are especially pronounced, since the 
scaling dimension xp{(3) turns out to have a vertical tangent at j3c- This might also 
contribute to the relatively poor outcome of the FSS analysis of the peak heights of 
the polarizability. With respect to the values of the effective correction exponents 
LO found for (3 < 13c (cf. table , we note comparing to the critical point behaviour 
that the nature of the corrections seems to be rather different in both cases, such that 
the effective correction exponents and amplitudes exhibit fast variation as the critical 
phase is entered, which is again related to the singularity of xp{f3) at pc- 

Finally, from deliberately reducing our simulation data set, we note that including 
lattice sizes only up to, e.g., L = 128, most of the estimates for /3c, 7/1^, P/i^, xp{(3) and 
p are not found to be compatible with the asymptotic results in terms of the statistical 
errors. Thus, consideration of large system sizes is crucial here for the resolution of 
scaling corrections, see also reference |0. This explains troubles experienced in the 
numerical analysis of the F model on a particular, annealed ensemble of fluctuating 
quadrangulations, which due to their intrinsic fractality only allow simulations of 
lattices with rather small effective linear extents . For a more detailed investigation 
of the thermal scaling properties an analysis involving measurements of the finite-size 
correlation length would be valuable. This, as well as the examination of the critical 
phase below the free-fermion point /3j, is left to a future investigation. 
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